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We present a novel treatment for calculating the in-medium quark condensates. The advantage 
of this approach is that one does not need to make further assumptions on the derivatives of model 
parameters with respect to the quark current mass. The normally accepted model-independent 
result in nuclear matter is naturally reproduced. The change of the quark condensate induced by 
interactions depends on the incompressibility of nuclear matter. When it is greater than 260 MeV, 
the density at which the condensate vanishes is higher than that from the linear extrapolation. For 
the chiral condensate in quark matter, a similar model-independent linear behavior is found at lower 
densities, which means that the decreasing speed of the condensate in quark matter is merely half of 
that in nuclear matter if the pion-nucleon sigma commutator is six times the average current mass 
of u and d quarks. The modification due to QCD-like interactions is found to slow the decreasing 
speed of the condensate, compared with the linear extrapolation. 
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I. INTRODUCTION 

The spontaneously breaking of chiral symmetry in 
Quantum Chromo dynamics (QCD) has been an interest- 
ing topic for nuclear physics It is expected that this 
broken symmetry would be restored in a medium, and 
several signals have recently been suggested for its mani- 
festation in nuclear matter, such as the spectral enhance- 
ment near the 2Mt^ (M^r is the pion mass) threshold 0, 
the systematic appearance of the parity doublets in the 
missing-states region Q , and the vector manifestation jjj 
etc. It is also shown that the partial restoration of chiral 
symmetry has strong effects on the nuclear saturation Q 
and the nuclear equation of state j^] . 

One of the most relevant quantity to the chiral symme- 
try restoration is the scalar quark condensate {qq). The 
zero value of the condensate can be strictly proven to be 
a necessary condition for exact chiral symmetry. But the 
sufficiency is not available presently. In fact, Birse has 
shown that the pieces arising from low-momentum vir- 
tual pions should not be associated with chiral restora- 
tion |2J. However, more recent investigations are in favor 
of the hypothesis that the quark condensate is a leading 
order parameter . 

The quark condensate in vacuum (qq)o is indeed a vast 
quantity (~ —10^ MeV'^). However, this value will be 
modified in a medium. The change produced by the 
deuteron has recently been investigated 0. There are 
also considerable investigations on the modification due 
to nuclear and/or quark matter [ul 13 113, IH, EE EE 
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ESElEllpElElEllS. For more references, see 

relevant review articles [l|, 123 ■ 

It has long been speculated that strange quark matter 
(SQM), rather than the normal nuclear matter, might 
be the true QCD ground state [l^, or absolutely stable 
in mm HE The critical density of SQM depends 
on the quark mass scaling 's^l, while the quark mass 
scaling is closely linked to the quark condensate in quark 
matter [sEl- Therefore, the study of in-medium quark 
condensates is meaningful independently of whether it 
can serve as an order parameter of chiral restoration. 

Originally, the interacting part of the quark mass was 
taken to be inversely proportional to density "29", "so", "sj 
Is^ . That parametrization has no derivation, and leads 
to inconsistencies in the thermodynamic properties of 
SQM j^. Later, a straightforward derivation method 
is provided, and the scaling is changed to be inversely 
proportional to the cubic root of the density [s^ . 

One of the most popularly used methods to calcu- 
late the in-medium quark condensate is the Feynman- 
Helmann theorem. The main difficulty in this formalism 
is that one has to make assumptions on the derivatives 
of model parameters with respect to the quark current 
mass. To bypass this difficulty, we will apply a similar 
idea as in the study of strange quark matter by defining 
an equivalent mass. A differential equation which deter- 
mines the equivalent mass will be derived in detail, then 
it is solved within a QCD-like inter-quark interaction. 
The normally accepted linear behavior of the condensate 
in nuclear matter is naturally reproduced and extended 
to quark matter. However, the decreasing speed of the 
condensate in quark matter is merely half of that in nu- 
clear matter if the pion-nuclcon sigma term is 6 times 
the average current mass of light quarks. At higher den- 
sities, the decreasing speed is usually slowed, compared 
with the linear extrapolation, due to interactions. 



2 



We organize the present paper as follows. In the sub- 
sequent Sec.^ we describe the equivalent mass method, 
and apply it to the investigation of the quark condensate 
in quark matter. An equation to determine the equiva- 
lent mass is derived, solved analytically at lower densi- 
ties, and solved numerically with a QCD-like inter-quark 
interaction at higher densities. To test the correctness of 
the equivalent mass method, we apply it to nuclear mat- 
ter in Sec. IIIII The model-independent linear behavior is 
reproduced and the higher order's contribution is studied 
in the mean-field level. Sec. II VI is a short summary. 



The basic idea of the mass-density-dependent model of 
quark matter is that the system energy can be expressed 
in the same form as a proper non-interacting system. The 
strong interaction between quarks is included within the 
appropriate variation of quark masses with density. In 
order not to confuse with other mass concepts, we refer 
to such a density-dependent mass as an equivalent mass. 
Therefore, if we use the equivalent mass m, the system 
Hamiltonian density should be replaced by an Hamilto- 
nian density of the form 



(5) 



II. THE CHIRAL CONDENSATE IN QUARK 
MATTER 

Applying the Feynman-Hellmann theorem, respec- 
tively, to the vacuum |0} and the state |nb} with baryon 
number density Tib, and then taking the difference, as 
has been done for nuclear matter j one can obtain the 
following expression for the relative chiral condensate in 
quark matter: 



{m)n 



1 



Sn* dniQ ' 



(1) 



where e is the energy density of quark matter above the 
vacuum, mo = 7.5 MeV is the average current mass of u 
and d quarks, n* is a constant related to the pion mass 
= 140 MeV and the pion decay constant i^^ — 93.2 
MeV by 



3 dmo 



0.98 fm" 



(2) 



Here we have used the GellMan-Oakes-Renner rela- 
tion m 



(3) 



and the notations (0)nb = ("-blOlnb) and {0)o = (0|O|0) 
for an arbitrary operator O. 

From Eq. (0), the usual procedure to calculate the 
in-medium quark condensate will be to find an expres- 
sion for the energy density, carrying then the derivative 
with respect to the current quark mass. However, it is 
usually difficult to know the current quark mass depen- 
dence of model parameters, e.g. meson masses, coupling 
constants, and cutoffs etc., without full understanding 
how effective models are related to the underlying theory 
(QCD). To bypass this difficulty, we suggest an equiva- 
lent mass method bellow. 

The QCD Hamiltonian density for the two- flavor case 
can be schematically written as 



Hqcd ^ Hi, + 2moqq + Hi, 



(4) 



where H\. is the kinetic term. Hi is the interacting part, 
and mo is the average current mass of u and d quarks. 
Terms breaking flavor symmetry have been ignored. 



where m is the equivalent mass to be determined. Obvi- 
ously, we must require that the two Hamiltonian densities 
iJoqv and Hqcd have the same expectation value for any 
state 1^), i.e.. 



(6) 



Applying this equality to the state jrib) with baryon num- 
ber density n^, and to the vacuum state |0), respectively, 
taking then the difference, one has 



(-ffcqv)nb ^ {Hcqv)o — (-ffQCD)nb " {Hqcd)o- 



(7) 



Since we are considering a system with uniformly dis- 
tributed particles, or in other words, the density nb has 
nothing to do with the space coordinates, we can write 
(vE'|m(nb)gq|5') = m(nb)(^'|gg|5'). This equality is espe- 
cially obvious if we consider it in terms of quantum me- 
chanics: IvE") is a wave function with arguments rib and 
coordinates, the expectation value is nothing but an in- 
tegration with respect to the coordinates. Therefore, if 
Ub does not depend on coordinates, the function m(nb) 
is also a coordinate-independent c-number, and can nat- 
urally be taken out of the integration. However, if rib is 
local, the case becomes much more complicated and we 
will not consider it here. 

Now we can solve Eq. Q for m, and accordingly obtain 



mo 



{Hl)rn, - {Hi)o _ 



— mo + mi. 



(8) 



Therefore, considering our quarks as a free system, i.e. 
without interactions, while keeping the system energy un- 
changed, quarks should acquire an equivalent mass of the 
form shown in Eq. ||SJ). From this equation we see that 
the equivalent mass m includes two parts: one is the 
original mass or current mass mo, the other one is the 
interacting part mi. Obviously the equivalent mass is a 
function of both the quark current mass and the density. 
At finite temperature, it depends also on the tempera- 
ture as well. Here we consider only zero temperature. 
Due to the quark confinement and asymptotic freedom, 
one may naturally expect 



lim mi 



oo 



and 



lim 

rib — *C)C 



mj — 0. 



(9) 
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Because the Hamiltonian density H^qv has the same 
form as that of a system of free particles with equiva- 
lent mass m, the energy density of quark matter can be 
expressed as 



9 
27r2 



^ \//c2 + k^dk = Smn^F , (10) 



where g — 2(flavor)x3(colors)x2(spins) = 12 is the de- 
generacy factor, 



kl = — TT^TT-b 
9 



1/3 



(11) 



is the Fermi momentum of quarks. 

Substituting Eq. H1(J|I into Eq. Q and carrying out the 
corresponding derivative, we have 



{qq)Q n* \mj drriQ 

where the chain relation 

de de dm 



(12) 



dmo dm, dm^ 

has been used. The functions F in Eq. (|10|l and / in Eq. 
((T^ are, respectively, defined to be 



F{x) EE - x^/x"^ + 1 (2x2 + 1) _ arcsinh(2;) /a;^,13) 
8 L J 

f{x) = ^ x^Jl!F^^ -\n[x+ \Jx'^ + l) jx^ (14) 
which have the following properties: 

F(x)-x^ = /(x), 



lim F{x) = 1, lim F{x) 



\im f{x) = 1, lim fix) ^. 



(15) 
(16) 

(17) 



Defining (i/i) to be the interacting energy density 
{Hi)rn, — {F[i)q and considering Eq. (0), the interact- 
ing part of the equivalent mass, mj in Eq. ((HJ, can be 
re- written as 



mi 



(Hi) /{in*) 



(18) 



1 - {m)n^l{m)Q 

Solving for the ratio {qq)n^l {(iq)o, this equation leads to 



{qq)n, ^ ^ 1 m 

(qq)o 3n* mi 



(19) 



According to Eq. the total energy density of the 
quark matter system can be expressed as 



e = A / Vk^+mok^dk + (Hi). (20) 
27r^ ' 

The first term is the energy density without interactions, 
the second term (Hi) is the interacting part. On the 



other hand, we have already expressed e in Eq. (|10() . So, 
replacing the e here with the right hand side of Eq. H10|l . 
then dividing by 3nb, we naturally get 



SUh \m J \ mo 

Substituting this into Eq. H19|) gives 



(qq)o n* m,i 



, fcf \ f kl 

mF I — — m^F 

m J \ mo 



(21) 



(22) 



Equating the right hand side of this equation with that 
of Eq. H12|l . one can easily get a fisrt order differential 
equation 



dm mF (kf/m) — moF (kf/mo) 



dmQ 



(m - mo)/ (fcf/m) 



(23) 



satisfied by the equivalent mass. Eqs. (|21|I - H23|I seem 
similar to the Eqs. (23)-(25) in Ref. jSg. They are in fact 
different by a factor of iVf, the number of quark flavors. 
This difference will lead to different numerical results. 

Now we have another formalism which can be used to 
calculate the in-medium chiral condensate. In principle, 
one can calculate the chiral condensate at any finite den- 
sity from Eqs. if^ with the aid of Eqs. ipSjl and 
if one knows some information on the interacting energy 
density (Hi) from a realistic quark model. The obvious 
advantage of this formalism is that one does not need 
to make assumptions on the current mass derivatives of 
model parameters. In the following, let's consider a sim- 
ple model as an example. 

If we denote the average distance between quarks by 
f, the interaction between quarks by v(mo,7ib)7 and as- 
sume each quark can only interact strongly with another 
A^o nearest quarks at any moment, because of the sat- 
uration property of strong interactions, we can link the 
interacting energy density (^^i) to density by 



(Hi) = -iVonbv(mo,r) 



(24) 



This is obtained as follows. Suppose the total particle 
number is N. The interacting energy due to particle i is 
iVov(mo,nb), accordingly the total interacting energy is 
^ -^ov(mo, rib) = ^NNov{mo, rib) (the extra factor 
1/2 is for the correction of double counting). Dividing 
this by the volume, we then have the interacting energy 
density in Eq. (|24|l where the average inter-quark dis- 
tance r is linked to density through 



1/3- 



(25) 



The average volume occupied by one particle is 
l/(3rib). If this volume is considered as a spherical ball 
of diameter f, one finds ^ = (2/7r)^/'^. If, however, this 
volume is considered as a cubic box of side f, one has 
^ = 1/3^/^. In the former case, one divides the system 
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into sub balls and place a particle at the center of each 
small ball. In the latter case, the system is divided into 
sub cubic boxes. Obviously there are unoccupied spaces 
between balls in the former case, so we take the latter 
value, i.e., C = 1/3^/^ 

Eq. (|24|l may not be an absolutely exact expression for 
the interacting energy density. Some other components, 
for example, the possible three-body interactions, have 
not been included. However, they should be of higher 
orders in density, and will not be considered here. To 
compensate for these ignorance, to some extent, one may 
regard the Nq as a free parameter. However, we will take 
No = 2 due to the fact that the quark has a trend to 
interact strongly with other two quarks to form a baryon. 
It will soon be shown that the relative chiral condensate 
is independent of the concrete A^o value at lower densities. 

Substituting Eq. into Eq. (|5T]) gives 



Nn 



mF 



moF 



Too 



(26) 



Keeping in mind the second equalities in Eq. jnj and 
in Eq. (|16|1 . we have 



lim v(too,?') 

nb— >oo 



0. 



(27) 



Note that we have said nothing about the form of the 
inter-quark interaction v(too, rib). 

At the same time, Eq. (|19|l should be compatible with 
Eq. H12|l . Comparing the right hand side of these two 
equations, we get 



, kf\ dm 

f \ ~] "^17^ 

ml omo 



No , 



(28) 



Because Eq. (I23II is a first order differential equation, 
we need an initial condition at toq = TOq to get a definite 
solution. Let's suppose it to be 



m{m*o,n\,) = m(nb). 
Usually, we will have 

v(TOo,nb)|mo=mS = v{f), 



(29) 



(30) 



where v{f) is the inter-quark interaction for the special 
value TOg of the quark current mass toq. 

Eq. H23|) is difficult to solve analytically. However, this 
can be done at lower densities. For convenience, let's first 
prove that the equivalent mass can be explicitly expressed 
through the inter-quark interaction as 



TO — mo + 



v(TOo,r) 



2 r f kt 

No J y ma 



c?V(mo,r) 
dnio 



(31) 



In fact, taking the derivative with respect to mo at 
both sides of Eq. gives 



dm 



./(M+^^lip^. (32) 
TO / omo \ mo I I omo 



Here Eq. H15() has been applied. Then substituting Eq. 
(I22J) into Eq. we immediately get Eq. 

At lower densities, the Fermi momentum fcf is small, 
so the function F{x) approaches to 1. Accordingly, Eq. 
(ESJ leads to 



Too 



No , 

— v(mo,r). 



(33) 



Replacing the left hand side of Eq. with this expres- 
sion, we get 



No dv 
2 dmo 



f 



\mo J 



= 1. 



(34) 



Integrating this equation under the initial condition given 
in Eq. (|30|l . we have 



v(TOo,r) = v{r) + / — 



mo 2 

No 



'-■'(-) 

\moJ 



dmo- (35) 



Then the interacting part of the equivalent mass is 
No 



TOl(TOo,nb) 



-v{r) 



1-/ 



Too 



dmo- (36) 



It is interesting to note that Eq. can, because of 
the dominant linear confining interaction, lead to a quark 

1/3 

mass scaling of the form toi oc 1 / which is consistent 
with the result in Ref. [33l |. 

Substituting the above results for v and toj into Eq. 
(fH^ gives 



(gg)«b 

{qq)o 



1 _ !^ 



(37) 



This result is very similar to the famous model- 
independent result for the quark condensate in nuclear 
matter 



{qq)c 



— with n* 

1* 



MlFl_ 
ctn 



(38) 



which was first proposed by Drukarev et al. (IL] . and 
later re-justified by many authors 0, 0, Q|. We will 
also reproduce it in the next section within the present 
approach. 

Both Eqs. H37|l and (|38|) mean a linear decreasing of 
the condensate at lower densities. The ratio of their de- 
creasing speeds is 



1/n* 



EL 

n* 



3too 
ctn 



(39) 



If one takes toq = (muo + mdo)/2 = (5 + 10)/2 = 7.5 MeV 
p3t and ctn = 45 MeV [stUMIsqII. i.e., the sigma com- 
mutator is about six times the average current mass, the 
decreasing speed in quark matter is merely half of that 
in nuclear matter. If, however, one takes for ctn a value 
of about three times TOq, as previously determined |40j |. 
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the two decreasing speeds are equal. More recent inves- 
tigations are in support of a bigger value for ctn 0, ■ 
In general, an explicit analytical solution for the con- 
densate is not available, and we have to perform numeri- 
cal calculations. For a given inter-quark interaction w(f), 
we can first solve Eq. (|26|l to obtain the initial condition 
in Eq. H29|l for the equivalent mass, then solve the differ- 
ential Eq. (|23|l , and finally calculate the quark condensate 
through Eq. 

There are various expressions for v{f) in literature, e.g ., 
the Cornell potential [i^, the Richardson potential p5l| . 
the so-called QCD potential s liol l47l | . and purely phe- 
nomenological potentials p8Ll49| |. etc. A common feature 
among them is that they are all flavor-independent. This 
independence is supported in a model-independent way 
by applying the inverse scattering approach to extract a 
potential from the measured spectra [s^. Let's take a 
QCD-like interaction of the form 



v{r) 



4a,(r) 



(40) 



The linear term af is the long-range confining part. It is 
consistent with modern lattice simulations and string in- 
vestigations |5ll Is^ . The second term incorporates per- 
turbative effects. To second order in perturbation theory, 
one has EE IH 



as{r) 



An 



6oA(r) 



1 - 



bi InA(f) 



A(r-) A(r-) 



where 535 



A(f)EEln[(fA^)-2+6] 



(41) 



(42) 



- 27Vf)/3, bi = [347V2 _ Arj(i37v2 _ 
(317Vc - 10A^f)/(9&o) for SU(iVc) and 



with 6o = (llA^c 
3)/A^c]/3, and 63 = 
Nf flavors. 

Besides these constants, there are three parameters, 
i.e., a, Ajns, and b. The QCD scale parameter is usually 
taken to be — 300 MeV. The value for the string 
tension a from potential models varies in the range of 
0.18—0.22 GeV2 [H, and we here take a = 0.2 GeV^. 
As for the parameter 6, we take three values, i.e., 10, 
20, and 30, in the reasonable range [H^. Because these 
parameters are determined from heavy quark experimen- 



tal data, the initial value 



is taken to be 1500 MeV 



which is compatible with the mass of heavy quarks. If 
we change this value to a smaller one, the global features 
of the condensate will not be changed significantly. Our 
numerical results are plotted in Fig. ^ and Fig.|21 

Figure n shows the density behavior of the chiral con- 
densate in quark matter. The straight line is the linear 
extrapolation of Eq. (|37|l . It does not depend on the form 
of the inter-quark interaction v{f), and so, at this mean- 
ing, 'model-independent'. The other three lines are for 
mo = 7.5 MeV, but for different b values, as indicated 
in the legend. At lower densities, the chiral condensate 
decreases linearly with increasing densities. When the 
density becomes higher, the decreasing speed is slowed. 
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FIG. 1: Density dependence of the quark condensate in quark 
matter. The straight line is the linear extrapolation of Eq. 
13711 . The other three lines are for mo = 7.5 MeV, but for 
different b values, as indicated in the legend. 



Figure [5] shows the current mass dependence of the 
chiral condensate in quark matter. From Eqs. H22(l and 
()23|l . it is not difficult to obtain 



d ( {qq). 



dmo V {qq)o 



71b 



F' 



[ dn* n* F* 

(43) 

[mF {{ki/m)) - moF {{kf/m,o))] /mi, 
and /g = fikf/mo). The first term in 



with F*(nh, mi) 
/* EE f{kf/m), 
the brace is due to the mo dependence of the vacuum 
condensate while the second term is purely due to the 
density. Here the numerical calculation only considers 
the second term to concentrate on the density effect. For 
the model Eq. H24|l with Eq. (|40|l . the condensate drops 
when the current mass becomes large. This effect is espe- 
cially obvious at higher densities. At lower densities, the 
condensate changes little with increasing current mass. 
This is understandable: because we have in fact assumed 
that the vacuum condensate or the quantity n* in Eq. 
is constant, the lower density condensate is nearly only a 
function of the density. If, however, the first term is neg- 
ative, i.e., the absolute value of the vacuum condensate 
decreases with increasing mo, as indicated in Ref. l56l|, 
the in-medium condensate will drops more rapidly. On 
the other hand, if the first term is positive enough, the 
in- medium condensate will go up with increasing mo. 

Finally in this section, we would like to comment that 
one can adopt any realistic model to calculate the den- 
sity dependence, and also current mass dependence if 
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section, we define an equivalent Hamiltonian density 
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FIG. 2: Current mass dependence of the quark condensate 
in quark matter. 



necessary, of the in-medium quark condensate without 
making any further assumption by following the formal- 
ism presented here, although we adopted a simple QCD- 
like interaction as an example. The procedure is as fol- 
lows. One first has the interacting energy density (Hi) 
from the selected realistic model for a special current 
mass value, gets the corresponding initial equivalent mass 
771 (nb) in Eq. H29|l from Eq. H21(l . obtains the equivalent 
mass 771(777.0, rib) by solving Eq. (j^ . calculates then the 
quark condensate through Eq. 122|l . 



III. THE QUARK CONDENSATE IN NUCLEAR 
MATTER 



In order to examine the correctness of the equivalent 
mass method, we now apply it to the investigation of the 
quark condensate in nuclear matter. 

To make a transit from quark matter to nuclear mat- 
ter, we suppose that nuclear matter can be well de- 
scribed by some effective Hamiltonian density which can 
be schematically expressed in the form 



(44) 



where Mn = 938.926 is the mass of nucleons in free space, 
Tik is the kinetic term, Tii is the interaction term between 
nucleons. We also ignore the nucleon mass difference. 
In principle, the following equality should hold for any 
nuclear state \p) with nucleon number density p: 



(pIHnmIp) = {p\Hqcb\p)- 



(45) 



Here Hqcb represents the QCD Hamiltonian density in 
Eq. Proceeding in the similar way as in the preceding 



Tic 



T^k + 2MNN, 



(46) 



where M is an equivalent mass to be determined by the 
requirement that the energy of this Hamiltonian is equal 
to the original interacting system, i.e. 



(plHcqvlp) = (pIHnmIp)- 



(47) 



Applying this equality to the nuclear state \p) and the 
vacuum state |0), respectively, and then taking the dif- 
ference, we have 

(p|Heqv|p) - (0|Heqv|0) = (pIT^nmIp) " (01Hnm|0). (48) 
Solving for M, this equation gives 

2 [(p|iViV|p} - (0|iViV|0)] ^ ' 
= M^ + Mi. (50) 

Obviously, Mi comes from the interaction between nu- 
cleons. It is generally a function of both the density p 
and the quark current mass rrio. If the density p is not 
too high. Ml is empirically much smaller than Mn- For 
example, the binding energy per nucleon is less than 2% 
of the free nucleon mass at the nuclear saturation density. 
Because no interaction exists when the density becomes 
zero, we naturally expect 



lim Mi = or lim M = M^. 



(51) 



Using the equivalent mass, the energy density of the 
nuclear matter can be expressed as 



27r2 



+ M2 p'^dp, 



(52) 



where = 2(isospins)x 2 (spins) = 4 is the degeneracy 

1/3 

factor, pf ~ (Gtt'^ p / gj^'j is the Fermi momentum of 
nucleons. 

From the application of the Hellmann-Feynman the- 
orem and the GellMann-Oakes-Renner relation, as was 
done in Ref . (l^j , one can easily get 



(gg)p 



1 deN 
377* dmo 



(53) 



Here en is the energy density of nuclear matter, and the 
quantity 77* is given in Eq. 

Substituting Eq. ifS^ into Eq. (|^ . we have 



{qq)o 



1 d£N dM 
377* dM dmo ' 



3n*-' \mJ 



d 



dmo 



(54) 

[A/n + Mi] . (55) 



p*-' \M 



mo dMi 



(56) 
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0.2 0.3 0.4 

DENSITY (fm ') 



FIG. 3: Density dependence of the equivalent mass M and 
its interaction part Mi in the Dirac-Brueckner approach. It 
is obviously shown that Mi is much smaller than Mn. 



In obtaining the third equality, we have used the defini- 
tion of the pion-nuclcon sigma term [ssj 



(Tn 



mo— or 



drriQ 



drriQ 



ctn 
Too' 



(57) 



The physical meaning of the equivalent mass is that 
nucleons should have the equivalent mass if the system 
is free but with unchanged energy. It is not difficult to 
prove that such an equivalent mass always exists. In 
principle, if one obtain the energy density with some 
realistic models or even from QCD in the future, the 
equivalent mass M can be obtained by solving the equa- 
tion H52() . For example. Fig. |21 shows the result from the 
Dirac-Brueckner method with Bonn-B potential 

It should be emphasized that the equivalent mass M 
is different from the conventional effective mass M* of 
nucleons. For example, at lower densities in Quantum 
Hadrondynamics (QHD) the equivalent mass and 

effective mass can be expressed as 



M* 
M 



Mn 



9 a 

Ml 



(58) 
(59) 



where Ma = 550 MeV and M^ = 783 MeV are, respec- 
tively, the masses for the a and w mesons, and ga = 10.2 
and gc^ = 12.55 are the corresponding coupling constants. 

It is obvious from Eqs. (|55)l and (|551) that the con- 
ventional effective mass includes merely the contribution 
from scalar mesons while the equivalent mass also in- 
cludes the contribution from other mesons. 



Equations H58|l and (|59() are obtained as such. The 
QHD Lagrangian is 

L = ^[-f''{id^~g^uj^)~{MN-ga<y)]'if + ld^<yd^a 



1 - .2 2 1 



(60) 



At the mean-field level, we have the following energy den- 
sity m 



eQHD 



■Pi 



9n 

lo 

{Mm - M*Mf + 



2M,^ 



(61) 



where the effective mass AI* is determined by the self- 
consistent equation 

M^2n^J, ^p2 + M*' ^ ' 

obtained by minimizing the energy density Eq. H61|l . 

Taking the approximation yjp + M*2 « M* in Eq. 
(|^ gives Eq. immediately. As for Eq. ifS^ , it should 
be solved from 



eN{M,p) = eQHD- 



(63) 



Here the eM{M,p) is the right hand side (r.h.s.) of Eq. 
(|52|l and egno is the r.h.s. of Eq. (|61|l . Taking the ap- 
proximation 



^/f~+AP « M and \J + M*^ w M* (64) 
in Eq. (|63|l . Eq. H59|) is accordingly obtained. 



1100 
1000 



900 --.^ 

800 - 

^ 700 - 

■J 600 - 

^ 500 - 

^ 400 - 

300 - 

^ 200 - 

100 - 

— ■ 

-100 L- 



M 

m' 

■M, 



0.0 



0.1 



0.2 0.3 
DENSITY (fm ') 



0.4 



0.5 



FIG. 4: The effective mass M* , the equivalent mass M, and 
the interaction mass Mi vary with density in QHD. 

Naturally, Eqs. H58|l and (|59|l are valid only to first 
order in density. For accurate relations, we have solved 
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Eqs. (|62|l and H63|) numerically. The results are plotted 
in Fig.H 

From Fig. 13 and Fig. 01 we see clearly that Mi is much 
smaller than M^. Up to three times the normal nuclear 
density, Mi is less than 6 percent of A/n- Also, when the 
density does not exceed the deconfinement density. Mi 
must be a continuous function of density. Therefore, we 
ignore the second term in the bracket of Eq. H55|l. and 
get 



(65) 



Because the second term in the bracket of Eq. H5t)|) has a 
factor mo/cTN which is much smaller than 1, we can also 
ignore the corresponding term and get the same expres- 
sion. 

Due to the smallness of the ratio pf/M at lower densi- 
ties, the function f{p{/M) approaches to unity. Accord- 
ingly, Eq. linni) becomes Eq. although Eq. (^3)) is 
in principle more accurate. Naturally, if one wants to 
study the higher density behavior of the quark conden- 
sate in nuclear matter, the function f{pf/M) should be 
kept. Also at higher densities, the omitted interacting 
part Ml should be taken into account. In the following, 
we try to consider it in the mean-field level. 

As has been found by Malheiro et al. 16], the quark 
condensate can also be linked to the trace of energy- 
momentum tensor of nuclear matter, i.e. 



= 1 - 



1 eN-3P 



(66) 



To continue our discussion, we now adopt this equality 
directly in our present model. Because it is justified by 
Malheiro et al. |l6j| in the mean-field approximation of 
Walecka, ZM, and ZM3 models, our treatment in the 
following are expected to be at least valid at mean-field 
level. In our case, the energy density eN is related to the 
equivalent mass M by Eq. 1)52(1 which, after performing 
the integration, gives 



{M,p)=MpF(^^) 



(67) 



Here the function F is the same as in Eq. (|13ll . 
The pressure P can be obtained from Eq. (|67() . i.e. 



Pf\ , 2 9M ^ f Pi 



where the function / is given by Eq. (|13|l . the function 
G is defined to be 

1 



G{x) = - x{2x^ - 3)Va?^+T-|- 3arcsinh(x) /x^. (69) 
8 L J 

The second term on the r.h.s. of Eq. H68|l is from the 
density dependence of the equivalent mass. In obtaining 
the pressure P, one should be very careful. Otherwise, 



inconsistencies will appear in this model. For thermody- 
namic details, see our recent publication ^2J|- 
Combining Eqs. H55|) and H66|l gives 



P_j. |^_P£^| dM _ 1 eN 



3P 



Sn*' \MJ dmo p* Mn 
substituting Eqs. and (EHl into Eq. {701), we have 



(70) 



dM 



M 
Mn 



F(pi/M) - 3G{pt/M) 3p dM 



3n* dmo 
Because of the equalities 
F{x) - 3G{x 



fiPf/M) 



M dp 



P 



Eq. 171|) can be cast to 



= 1 and J— = ^ 

3n* 



mo 

f7N ' 



mo dM 3p dM _ M 
CTN i9mo Mn dp Mn 



(71) 
(72) 

(73) 



This partial differential equation is linear and of first 
order. To find a special solution for our purpose, let's 
expand the equivalent mass M in a Taylor series at zero 
density, i.e. 



M{p, mo) = c^p' 



co + cip + C2p^ + C3p^. (74) 



In principle, this expansion can have an arbitrarily large 
number of higher order terms in density. But, for the 
moment, we only consider terms up to the third order in 
density. We will discuss this problem at the end of this 
section. 

Substituting Eq. into Eq. (O leads to 



mo 
(Tn 



Mt^ 



E 



lap 



(75) 



where c/ = dci/dmQ. Comparing the coefficients in Eq. 
(|7^ . we have 



ctn 3i — 1 
mo Mn 



-Ci, 



or explicitly. 



Co 
ci 

/ 

C3 



(IN Co 

mo A/n ' 

O'N Ci 

mo Mn' 

C2 

mo Mn' 

0"N C3 

mo Mn ' 



-2 
-5 



(76) 

(77) 
(78) 
(79) 
(80) 



From Eqs. l(5T|l and (fZ^ . we can easily get cq — Mn. 
In this case, Eq. (|77|l becomes cq' = af^/ruQ, which is 
consistent with Eq. (|57|l . 
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It should be pointed out that the coefficients Ci in the 
Taylor expansion l|74|) are mQ-depcndcnt. In fact, we can 
integrate Eq. (f7?)|l and get 



c.(7.5) 



exp 



7.5 Mn mo 



3j-l 



(81) 



where Ci(7.5) is the value for Ci at mp = 7.5 MeV which is 
the value we take for the current mass of u and d quarks. 
Now let's substitute Eq. into Eq. or (j^ and 

get 



(gg)p 
{qq)Q 



1 



1 



mo 
(Tn 



(82) 



On applying Eq. (|75|l . we finally have 

3 



(gg)p 



1 



(83) 



For figure-plotting convenience, let's change the coef- 
ficients Ci (i=l,2,3) to ai by defining 



-Ci or c, 



(84) 



A/n Po 
Accordingly, Eqs. ifT^ and become, respectively, 

# = l + a,^+aJP-)\aJP-)\ (85) 
A/n Po Vpo/ VPo/ 



p* 



1 - E(3^ - 



86) 



To obtain the values of the coefficients ai, 02, and 03, we 
can use the following known experimental properties of 
nuclear matter 



eN(M,p) 
P 



dp 

^pI£ 



P=Po 
eN(M,p) 



eN(A^,p) 
P 



A/N + £^b, 

= 0, 
= /ST. 



P=PO 



(87) 
(88) 
(89) 



Here po = 0.17 fm^"^ is the normal nuclear saturation 
density, = —16 MeV is the nuclear binding energy. 
As for the incompressibility it is not so definite, 
maybe ranging from 150 to 350 MeV. The expression 
eN(A/, p)/p = MF{pi/M) is the energy per nucleon. For 
a given K, we can solve Eqs. H87I89|I and get the values 
for A/(') (i==0,l,2): 



M(") EE M\ 



P=PO ' 

dM 

P=PO 



(90) 
(91) 

(92) 



0,06 
0,05 
0,04 
0,03 
0,02 
0,01 
0,00 

CO" 

-0,02 

^ -0,03 
-0,04 
-0,05 
-0,06 
-0,07 
-0,08 



o a, 

- <^-- ^3 

Pc 

Pd 



o 00-0 



r O O O O 



O O O 



O O O 



O O O 



250 

K (MeV) 



FIG. 5: The coefficients ai, a2, and 03 vary with the incom- 
pressibility K. The dotted line gives the corresponding crit- 
ical density pc at which the quark condensate becomes zero. 
The full line is the deconfinement density pd after which the 
deconfinement phase transition is expected to occur. 



Substituting the expansion (|85|l into the right hand side 
of Eqs. (|9()I92(I and then solving for 01,02, and a^, we 
have 



Ol = 



0-2 



-1 

Mn 
1 

A^ 

-1 

Mn 



3(A/n - A/(°)) + 2A/(^Vo - ^A/^^Vo^ 
3 (A/n - A/(°)) + 3A/(iVo - M^^Vo^ 



A/n - A/(°) + A/(i) 



Po 



A/(2)po2 



q3) 
,(94) 
(95) 



P=PO 



The quark condensate can then be calculated by Eq. H86|l . 

In Fig.|Sl we give, on the left axe, the coefficients ai, a2, 
and 03 as a function of K. Another two lines on the right 
axe are to be explained a little later. 

Fig.l^shows the equation of state of nuclear matter for 
typical incompressibilities. The open circle corresponds 
to the saturation which is a mechanically stable point, 
so every line passes it by. The points marked with a 
solid circle are also in mechanical equilibrium, but they 
are not stable. When the density is between the open 
and solid circle, the pressure is positive. However, when 
the density exceeds the full circle, the pressure is drasti- 
cally negative, so the system will automatically contract 
and the deconfinement phase is expected to appear. We 
therefore use pd to remember this density and plot it in 
Fig. |31 with a solid line on the right axe. We can also see 
that the incompressibility K should be no less than 265 
MeV because nuclear matter will be less bounded than 
16 MeV on the right side when K < 265, i.e. the differ- 
ence of the energy corresponding to the solid and open 
circles is smaller than 16 MeV if K is less than 265 MeV. 
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FIG. 6: Equation of state of nuclear matter for different in- 
compressibilities in the present approach. The open circle 
corresponds to the saturation. After the solid circle, nuclear 
matter will automatically contract, and the deconfinement 
phase is expected to occur. 



The quark condensate is given in Fig.|7|for typical val- 
ues of K . For K > 265, the critical density pc where the 
quark condensate becomes zero is generally bigger than 
that obtained by extrapolating the model-independent 
result of Eq. H38|l (indicated with the word 'Linear' in 
the figure legend). The variation of the critical density 
with the incompressibility is shown in Fig. [SI with a dot- 
ted line on the right axe. A noticeable feature is that, for 
K > 326, there appears a hill where the chiral symmetry 
breaking might be enhanced. 

As mentioned before, the Taylor expansion Eq. H74|l 
can have infinitely large number of terms in density. 
However, it has been obviously shown in Fig. |31 that 



|ai| > \a2\ > laal 



(96) 



This is an indication that the higher the term's order is, 
the less important its contribution. This is why we have 
only taken to the third order in the expansion Eq. (|74f) . 
However, it is very easy to extend to even higher orders. 
For this purpose, we should first express the equivalent 
mass as 



FIG. 7: The Quark condensate in nuclear matter to its vac- 
uum value in the present model. When the incompressibility 
K is in the range 265-326 MeV, the condensate decreases 
monotonously with increasing density. When K is bigger 
than 326 MeV, there appears a hill where the chiral sym- 
metry breaking might be enhanced. In both case, the critical 
density for the condensate becomes zero is higher than the 
linear extrapolation. 



where J is an integer. It's value depends on to which 
order we would like to expand. In the case we have just 
treated, J = 3. 

The initial values for 



dp' 



(z = 0,---, J-1) (98) 



are solved from 



(3po)' 



dp* 



eN(M,p) 



P 



p—rhoQ 



= K,, (i = 0.--J-l), 
(99) 

where Kq = Mt<, + E^, Ki ^ 0,K2 = K. If J > 3, we 
should also know the information on if 3, K4 et al.. The 
coefficients {i = 1 . . . J) are obtained by solving 



M{p, mo) = ^ 



CtP 



(97) 



E 



{i-ky. 



c,po^ = M^^\ {k^O---J -I). (100) 
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And the quark condensate is 




(101) 

Here the relation between and Ci is still given by Eq. 

IV. SUMMARY 

We have presented a new formalism for the calcula- 
tion of the in-medium chiral condensates. The key point 
of this scheme is to obtain an equivalent mass by solv- 
ing a differential equation. The main advantage is that 
one does not need to make further assumptions on the 
derivatives of model parameters with respect to the quark 
current mass. 

As an application of this method, we have successfully 
reproduced the normally accepted model-independent re- 
sult in nuclear matter. We also showed that there is a 
similar expression for the chiral condensate in quark mat- 
ter. It is found that the pion-nucleon sigma term is six 
times the average light quark current mass if the two 
decreasing speed is equal. 

For quark matter at higher densities, we adopted a 
QCD-like interaction to solve our differential equation, 
which shows that the condensate decreases linearly at 
lower densities. However, the decreasing speed is slowed 



at higher densities, compared with the linear extrapola- 
tion. For some values of the parameters, it can finally 
vanish, while for other values it may saturate and even 
increase. 

We have also shown that the incompressibility K 
should not be less than 265 MeV to give the correct 
nuclear binding energy in this model. When K is in 
the range of 265 to 326 MeV, the quark condensate 
monotonously decreases with increasing density. How- 
ever, the critical density where the condensate vanishes 
is, generally, higher than that obtained by linear extrap- 
olation. When K is greater than 326 MeV, a hill for the 
condensate as function of density will appear, enhancing 
the chiral symmetry breaking, vanishing then for even 
higher densities. 
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